MA/UC3M/05/97 
dc motion of ac driven sine-Gordon solitons 



Niurka R. Quintero^ and Angel Sanchez* 
Grupo Interdisciplinar de Sistemas Complicados, Departamento de Matemdticas, 
c/ Butarque 15, Universidad Carlos III de Madrid, E-28911 Leganes, Madrid, Spain 

(February 1, 2008) 



Abstract 



We investigate the possibility of dc soliton motion sustained by pure ac 
driving in the sine-Gordon model. We show by means of the collective co- 
ordinate formalism that ac driving induces a net dc velocity whose modulus 
and direction depend on the driving phase. Numerical simulations of the full 
sine-Gordon equation confirm the correctness and accuracy of this prediction. 
Non trivial cases when dc soliton motion is transformed into oscillatory as 
well as the effects of damping are analyzed. Our results settle a long stand- 
ing issue about the existence and characteristics of this phenomenon, whose 
possible appearance in other systems is also discussed. 

PACS number(s): 03.40.Kf, 74.50.+r, 85.25.Cp 
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I. INTRODUCTION 



It is well recognized that the study of nonlinear equations and their solutions is of great 
importance in many areas of physics. Among these, the sine-Gordon equation has attracted 
the interest from mathematicians and physicists due to not only its complete integrability 
but also its ubiquity as a model of nonlinear physical phenomena. Indeed, the sine-Gordon 
equation is known to be a canonical model for a wide variety of physical systems when 
topological solitons are present, such as motion of dislocations in crystals [|^, charge density 
waves [0], solitons in magnetically ordered systems 0, epitaxial growth of thin films [Q, 
fluxon dynamics in long Josephson junctions [^], or DNA promoter dynamics 

One still unsolved issue regarding the application of external forces on systems repre- 
sented by the sine-Gordon model regards the effects of pure ac driving and the possibility 
that it induces dc motion of solitons. To clarify this question, in this Communication we 
analyze the corresponding perturbed sine-Gordon equation 

(f)u-(f)xx + sm{(j)) = -(3(f)t + ef{t), (1) 

where subindices indicate derivative with respect to the corresponding variable, /30t is the 
usual damping term, and f{t) = sin{6t + 6o) is an external periodic force describing, for 
example, a long Josephson junction under the application of a uniform microwave field 
0. In the following, we will only consider ac driving (with or without dissipation) acting 
on a single soliton. We note that the influence of an ac force like the one we study here 
acting on breather solutions of the sine-Gordon equation was the subject of several papers 
by the end of the eighties, mostly from the view point of the use of a suitable frequency 
to stabilize breathers in the presence of dissipation P,^. The related problem of an ac 
driving in the presence of an additional dc driving has been already considered and it is 
quite well understood (see, e.g., |0,|lO|,|lT| and references therein), whereas the effects of pure 



dc driving were established in the seminal work of McLaughlin and Scott |T^. The issue of 
pure ac driving acting on solitons, however, has proven itself much more confusing, and a 
clear picture of soliton behavior under this kind of force was still lacking. 

To our knowledge, ac driven sine- Gordon solitons were first studied by Bonilla and Mal- 
omed [0 (and subsequently in for the Toda lattice with similar results), who claimed 



that ac driving could support dc motion of solitons in a discrete (in space) sine-Gordon 



model, dissipation being crucial for this phenomenon. Unfortunately, as was shown in [|15 
this result was incorrect because of two reasons: First, it was obtained by means of a neces- 
sary, but not sufficient, existence condition, and second, the authors used simultaneously and 
carelessly both the discrete and continuum limits of the model. Numerical simulations con- 



firmed that the dc motion predicted in |jT3[ did not take place . Recently, another study 
found a parameter region for the discrete sine-Gordon model (i.e., the Frenkel-Kontorova 
chain) where ac driving can induce dc motion. Such behavior can only occur in the discrete 
model because its two ingredients arise from discreteness: the Peierls-Nabarro barrier and its 



associated frequency [|T6| . Hence this is a much more complicated (indeed, chaotic) process 
involving attractor competition and has nothing to do with the proposal of [jl3|. 

In this work, we show that dc motion of solitons induced by pure ac driving in the 
continuum sine-Gordon equation is indeed possible: In the next Section, we use a collec- 
tive coordinate approach to compute the soliton motion, finding that its mean velocity is 
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a function of the phase of the driving. We verify our results by means of numerical sim- 
ulation, the comparison turning out to be satisfactory with an almost perfect agreement. 
The final Section contains a discussion of our results and comments on related subjects and 
generalizations to other systems. 



II. COLLECTIVE COORDINATE ANALYSIS 



In order to perturbatively study eq. (||) we will resort to the usual collective coordinate 
approach |ll7| , p!8| : If (3 and e are small, we can approximate the one-sohton solution of eq. 
(HI) by that of the unperturbed sine-Gordon equation (/3 = e = 0) with time dependent 
parameters, i.e., we will make the following ansatz: 



(x, t) = A arctan exp 



Xoit) = / u{t')dt' 
Jo 



X - Xo{t) - Xpit) 



(2) 



where the positive (negative) sign corresponds to a kink (anti-kink) solution. We note that 
X{t) = Xo{t) + Xo{t) [u{t)] has the meaning of the position (velocity) of the center of the 
soliton, hence the name "collective coordinate," and that the main assumption underlying 
this approximation is that radiation effects induced by the perturbation are neglected. This 
will be verified a posteriori by comparing with the numerical simulations. 



We now apply the method of McLaughlin and Scott [|T2| to obtain the equations of motion 
for the soliton center. The method, which they dubbed energetic analysis, is a rather simple 
one and amounts to compute the variation of the energy and momentum of the unperturbed 
sine-Gordon system due to the perturbation and the same quantity for a soliton solution, 
subsequently imposing compatibility. This is already a classic procedure and details can be 
found in [12,17,|l8|. The resulting equations of motion are (again, the ± sign corresponds to 



a kink or an anti-kink) 



du 



-i(l - u^){±n^l-v?tf{t) + A(3u), 



dX 



(3) 



We begin by discussing the dissipation-free case, i.e., we set /3 = in eqs. (|^). In this case, 
the first equation of (H) can be exactly solved, yielding 

u{t) = F[u{Q),5o\/{l + F[u{Q),5,ff'\ 

F[w(0), 5o] = , ± -r-Mos{5t + 5,) - cos(<5„)]. 

'I - m2(0) 45 
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± 



Equation (^) should now be integrated to obtain the final result, namely the soliton center 
motion. This cannot be done in general; however, when |'u(0)| << 1 and |||| << 1, a good 
approximation to the solution is given by 

X{t) = X(0) + [u{0) T ^ cos {5o)] t± 

Ad 

(5) 

7ie[sm{6t + 6o) — sin(5o)] 

From eq. (j^) it is evident that kinks and anti-kinks will move with net velocity {u{0) =F 
[it e cos {60)]/ 4:6} in a straight line, to which oscillations of frequency 6 are overimposed. 
Therefore, in this approximation, an ac driving of frequency 6 will almost always induce a 
dc motion of sine-Gordon solitons; indeed, only if 

^0) = ±!Li^ (6) 

solitons will remain oscillating around their initial position. It is important to note that this 
condition depends on the initial velocity of the soliton, u{0); if ■u(O) = 0, i.e., the soliton 
is initially at rest, the condition is simply 6q = (2n + l)vr/2, n = 0, ±1, . . . But even if the 
soliton is initially moving, a phase can be chosen such that this dc motion is stopped and 
transformed to an oscillation, a non-trivial and unexpected result. 

We note that the above predictions have been obtained for the case of a slowly moving 
soliton and a small driving strength. It is possible, however, to establish this result in general 
by means of an alternative method based on the study of the Hamiltonian. This calculation 



is somewhat more complicated, so we will report on in elsewhere [19| and here we will just 
quote the final result, namely condition for oscillatory motion becomes 

n(0) 7recos(5o) 

i.e., the only change with respect to the approximate condition previously derived is the 
(otherwise expected) appearance of the Lorentz factor 7 = y/l — u{Oy. 

To conclude our analytical calculations for this problem, we now turn to the dissipative 
case. Then, eq. (3) cannot be solved, not even for u. A slow motion approximation << 
1) yields 

u(t) = c exp(— /9t)=F 



T -^^^^-^[/3sm{5t + 60) -6 cos(5t + 5o)], 



(8) 



with 



evr 

c = u{0) ± ^^^^-^^[f3sm{So) -ScosiSo)]. (9) 

From the above, rather involved expression [or from the even more cumbersome one for 
X{t) which can be obtained by integrating (P)], the main conclusion we can draw is that 
(within the small velocity approximation) solitons will never exhibit dc motion except for a 
transient, after which they will reach a final oscillatory state around a point depending on 
the initial conditions. 
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III. NUMERICAL RESULTS 



All the results and predictions in the previous section involve some approximation, be- 
ginning with the collective coordinate hypothesis of negligible radiation. Therefore, they 
are meaningless unless verified by numerical simulations of the full perturbed sine-Gordon 
problem (|l|). Simulations were carried out by means of a standard fourth order Runge-Kutta 
algorithm with initial conditions given by an unperturbed sine-Gordon soliton, at rest 
or with velocity m(0), and boundary conditions (px{L = ±50, t) = 0. The conclusion of our 
numerical simulation program was that the predictions of the collective coordinate theory 
are very well verified by the full perturbed sine-Gordon equation. An example of the accu- 
racy of our analytical results is shown in Fig. |I], where we plot the results of simulations in 
the absence of dissipation for different values of the phase 60 and a soliton (actually, a kink) 
initially at rest. The high degree of agreement between theory and simulations is apparent, 
confirming that dc motion of ac driven kinks is indeed the usual behavior, except for very 
particular choices of the phase of the driving. 

The rest of the predictions are equally correct. Thus, Fig. ^ collects the outcome of 
simulations intended to verify the accuracy of the prediction that solitons with specific 
velocities can be stopped by a suitable choice of the driving phase. To verify this result, 
we proceed in the opposite way, i.e., we compute the critical velocity for a given value of 
the phase, and then we obtain the same value numerically. As can be seen from Fig. |^, 
the numerical value turns out to be u{0) = 0.1335, to be compared to the predicted m(0) = 
0.1347, i.e., the accuracy is better than 1%. Finally, Fig. 3 is an example of the behavior of 
the system when dissipation is present. We find numerically that, indeed, oscillatory motion 
is the only outcome of simulations for whatever value (not even small ones) of the initial 
velocity (which, in turn, agrees with the results in [p!5|] ). 



IV. CONCLUSIONS 

In summary, in this Communication we have analytically and numerically demonstrated 
the possibility of soliton dc motion induced by pure ac driving in sine-Gordon systems. 
The direction and modulus of the corresponding velocity have been shown to depend on 
the driving phase, which for a specific value (depending on the initial data) leads to pure 
ac motion even if the soliton was moving. The analytical results, obtained by means of a 
collective coordinate approach, agree with the numerical ones within an accuracy better than 
1%. We have also found that this effect is only possible in the absence of dissipation, and that 
the final steady state for all damping values is the oscillatory one. This conclusion finally 
clarifies the problem of ac driving effects on solitons and opens the way to its applications 
in different technological contexts, such as, e.g., Josephson junctions @J^. 

To close this section, some comments are in order about the existence of this phenomenon 
in other systems. We believe that the features we have found here are generic in the collective 
coordinate sense, i.e., for systems whose solitons exhibit a clear and robust particle-like 



behavior (see, e.g., ||17| , p!8| and references therein). However, care must be taken when dealing 
with solitons possesing inner degrees of freedom, such as the sine-Gordon breather or the 
0^ solitary wave, as the frequencies pertaining to these degrees of freedom will necessarily 
interact in a complicated manner with the driving frequency. It is clear that this competition 
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can give rise to chaotic phenomena, similar to that found in |T^, overruhng the simple 
scenario we have depicted here. On the other hand, if those frequencies are very different, 
our results should once again hold. Preliminary analytical and numerical investigations 
support this conclusion [119 
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FIG. 1. Verification of collective coordinate predictions in the absence of dissipation. Sim- 
ulation starts from a static kink located at X{0) = —5, and subject to an ac force given by 
— 0.02 sin(0.1t + 6q). (a) Sq = (upper curve), 5o = x (lo'^sr curve); notice that the direction of 
motion is opposite in both cases, (b) 6o = ^, critical value exhibiting oscillatory motion. Both in 
(a) and (b), solid line corresponds to numerical integration of eq. (jl]), dashed line is the analytical 
prediction (^). 
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FIG. 2. Stopping of a kink initially moving with velocity u(0) = 0.1335 (solid line) compared 
to the evolution of kinks moving with u(0) = 0.135 (upper dashed curve) and n(0) = 0.132 
(lower dashed curve). The ac driving had 5q = vr/G; other parameters are 5 = 0.1 and e = 0.02. 
The theoretically predicted critical velocity was m(0) = 0.1347, and the corresponding theoretical 
evolution lies on top of the numerical curve for u{0) = 0.1335. 
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FIG. 3. Verification of collective coordinate predictions with dissipation. Shown are two exam- 
ples of soliton motion subject to an ac force given by esin(0.1t + 5q), with e = 0.01 (upper curve) 
and e = 0.05 (lower curve); other parameters are P = 0.01, 6o = 0.1 and u{0) = 0.2. Solid lines 
correspond to numerical integration of eq. (||), dashed lines are the analytical predictions obtained 
by integrating (S). 
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